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A modeling of the normal state of the p-doped high-temperature superconductors (HTSC's) is 
presented. This is achieved starting from a more conventional metallic phase for optimal- and 
overdoping and passing via the underdoped to the insulating state by consecutive orbital selective 
compressibility-incompressibility transitions in terms of sum rules for the charge response. The 
modeling is substantiated by corresponding phonon calculations. Extending investigations of the 
full dispersion and in particular of the strongly doping dependent anomalous phonon modes in 
LaCuO, which so far underpin our treatment of the density response of the electrons in the p- 
doped HTSC's, gives additional support for the modeling of the electronic state, compares well with 
recent experimental data and predicts the dispersion for the overdoped regime. Moreover, phonon 
densities of states have been calculated and compared for the insulating, underdoped, optimally 
doped and overdoped state of LaCuO. From our modeling of the normal state a consistent picture 
of the superconducting phase also can be extracted qualitatively pointing in the underdoped regime 
to a phase ordering transition. On the other hand, the modeling of the optimal and overdoped state 
is consistent with a quasi-particle picture with a well defined Fermi surface. Thus, in the latter case 
a Fermi surface instability with an evolution of pairs of well defined quasiparticles is possible and 
can lead to a BCS-type ordering. So, it is tempting to speculate that optimal Tc in the HTSC's 
marks a crossover region between these two forms of ordering. 
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I. INTRODUCTION 

In our microscopic theoretical description of the charge 
response, electron-phonon interaction (EPI) and lat- 
tice dynamics of the high-temperature superconductors 
(HTSC's) in the framework of the linear response ap- 
proach we focus on the specific features of their solid- 
state chemistry. This is different from the overwhelming 
number of attempts relying on simplified models with 
strongly reduced model Hamiltonians. Nevertheless, this 
seems necessary since the phenomena observed in the 
cuprates hardly even occur in the other known com- 
pounds. One important feature neglected for example 
in the popular models of Hubbard type is the proxim- 
ity of the HTSC's, even in the optimally doped case, 
to the layered ionic crystals where besides short-ranged 
correlations strong long-range Coulomb interaction and 
strong (nonlocal) EPI exists. In the absence of full ab 
initio methods for the HTSC's to tell us what the correct 
starting point would be, such insight is essential, besides 
the proximity of the materials to an antiferromagnetic 
state. Models with strongly reduced electronic degrees 
of freedom may be useful for the discussion of theoretical 
ideas in general rather than for detailed calculations of 
the material specifica of the cuprates. 

The strong ionic nature of the HTSC's is modeled 
in our approach by an ab-initio rigid ion model (RIM), 
used as a reference frame for the local ionic, rigid charge 
response and the EPI, respectively. The important 
and characteristic nonlocal, non-rigid contribution to the 



electronic density response and the EPI in the HTSC's 
is calculated in terms of microscopically well defined 
charge- and dipole fluctuations (CF's and DF's). The 
latter are excitable on the electronic shells of the ions in 
the crystal. This is, because most of the charge in the 
HTSC's contributing to bonding is concentrated on the 
latter. Additionally, covalent metallic features of bond- 
ing are approximatively taken into accoun1jii£. In the 
calculations a sufficiently broad orbital freedom (Cu3d, 
Cu4s and 02p), the full three-dimensional Coulomb in- 
teraction as well as the short-ranged local repulsions, in 
particular the important on-site repulsion mediated by 
the localized Cu3d orbitals, is considered quantitatively. 
All the resulting couplings appearing in the dynamical 
matrix and the EPI are microscopically well defined and 
can be calculated. The results obtained so far agree well 
with the experimental phonon dispersion. For a recent 
review, see Ref. Ql 

The model for the electronic polarizability matri>si*^ 
II Kre '(q) of the optimally-, and overdoped metallic state 
of the HTSC's, describing the kinetic part of the charge 
response, is obtained in a first step from the Kohn-Sham 
parameters (usually interpreted as bandstructure) within 
the local density approximation (LDA) of density func- 
tional theory (DFT). The indices k, k' denote the or- 
bital degrees of freedom in the elementary cell of the 
crystal, i.e. Cu3d, Cu4s and 02p. q is a wavevector 
in the first Brillouin zone. However, such a model, ap- 
plied to LaCuO overestimates the coupling along the c- 
axis and thus underestimates the real anisotropy in the 
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cuprates. Modifications to improve the charge response 
along the ionic c-axis have been investigated in Ref. 00- 
Given the nearly two-dimensional electronic structure of 
the cuprates and thus a very weak interlayer coupling, 
electron dynamics and phonon dynamics will be on the 
same time-scale in a small region around the c-axis. This 
needs a non-adiabatic treatment with dynamical screen- 
ing of the bare Coulomb interaction and phonon-plasmon 
mixing in the metallic state of the HTSC'sA^. 

In the absence of a rigorous quantitative description 
for the insulating (undoped) state of the HTSC's , the 
model for H KK > has been designed to be consistent with 
the characteristic electronic charge response for an insu- 
lator with a gap in the charge excitation spectrum. This 
is achieved by fulfilling rigorous sum rules for the density 
response in terms of the polarizability matrix in the long- 
wavelength-limi1j2i£. These sum rules, see Section lTTTI can 
be considered as orbital resolved closed forms to represent 
the change of the charge response in a metal-insulator 
transition in terms of the electronic polarizability or the 
compressibility of the electronic system, respectively, as 
a primary tool to characterize the ground state. 

Finally, a suitable model is proposed for the under- 
doped state of the p-doped HTSC's, where a description 
of the low-energy excitations remains a theoretical chal- 
lenge and where the question how the metallic charge 
conduction emerges from the doped insulator is poorly 
understood. This model describes a novel metallic state 
with an insulator-like, incompressible charge response re- 
lated to the localized Cu3d orbitals at the Cu sublat- 
tice interspersed with a metallic, compressible charge re- 
sponse related to the more delocalized 02p orbitals at 
the oxygen sublattices of the CuO plane. In this way a 
partial ordering of the conducting carriers in the CuO 
plane is obtained. With regard to our modeling of the 
metallic state for the optimally and overdoped HTSC's 
such an ergodicity breaking delocalization-localization 
transition in the underdoped state in terms of the com- 
pressibility means a compartmentalization of configura- 
tion space. This means that some parts of configura- 
tion space cannot be approached or are hardly accessible 
to the particles. On the other hand, this should be ac- 
companied by a reciprocal compartmentalization of mo- 
mentum space, according to the uncertainty relationship 
Ax • Afc ~ 1. Indeed, angle-resolved photoemission spec- 
troscopy (ARPES) studies in underdoped cuprate a 6 i 7 i 8 , 
reporting a "two component" electronic structure with a 
pseudogap around (n, 0) and only part of the Fermi sur- 
face surviving as an "arc" around the node, point to com- 
partmentalization of momentum space too. The com- 
partmentalization of momentum space counteracts the 
reduction of the electron scattering rate brought about 
by the Pauli-principle in a Fermi-liquid, at least partially. 
Thus, the scattering rate and the life time of the carriers 
no longer will follow the well known scaling with energy 
and temperature of a Fermi-liquid. Consequently, the 
normal state should display non-Fermi-liquid properties 
in our model for the underdoped state of the cuprates, 



also seen in the experiment. In the corresponding su- 
perconducting state an increase of the kinetic energy as 
in BCS by the "smearing" of the Fermi surface would 
not be present because the latter is not developed. On 
the other hand, phase coherence in a phase ordering sce- 
nario would enhance the mobility of pairs and decrease 
the kinetic energy. Such a description of a delocalization- 
localization crossover of the electronic structure (e.g. due 
to doping) accompanied by a pseudogap phenomenon 
in terms of consecutive orbital-selective compressibility- 
incompressibility transitions in the presence strong non- 
local EPI may also be useful to model the charge correla- 
tions in other systems besides the cuprates characterized 
by a mixture of localized and delocalized orbitals. 

The article is organized as follows. In Section [H] the 
theory and modeling is briefly reviewed to set the frame. 
Section HTT1 presents our calculated results of the phonon 
dispersion for the underdoped and overdoped state of 
LaCuO. Moreover, the phonon anomalies related to the 
high-frequency oxygen bond-stretching modes and the 
phonon density of states for the insulating, underdoped, 
optimally doped and overdoped state of LaCuO are cal- 
culated. A discussion of the modeling of the electronic 
state in the cuprates consistent with the corresponding 
phonon calculations is provided. Finally, a summary will 
be presented in Section llVl 



II. SKETCH OF THE THEORY AND 
MODELING 

From a general point of view our treatment of the elec- 
tronic density response and lattice dynamics in terms 
of DF's and CF's can be considered as a microscopic 
(semi-ab-initio) implementation of the phenomenological 
dipole-shell model or the charge-fluctuation models, re- 
spectively. For a general formulation of phenomenolog- 
ical models for lattice dynamics that use localized elec- 
tronic variables as adiabatic degrees of freedom, see for 
example Ref. 0. This formulation covers shell mod- 
els, bond-charge models and charge-fluctuation mod- 
els. While in this approach the coupling coefficients 
are treated as empirical fitting parameters, the essential 
point in our microscopic scheme is that all the couplings 
can be calculated. 

In the following a survey of the theory and modeling is 
presented. A detailed description can be found in Ref. [B| 
and in particular in Ref. ^3 where the calculation of the 
coupling parameters is presented. 

The local part of the electronic charge response and 
the EPI is approximated in the spirit of the quasi-ion 
approachii by an ab initio rigid ion model (RIM) taking 
into account covalent ion softening in terms of (static) ef- 
fective ionic charges calculated from a tight-binding ana- 
lysis (TBA). In addition, scaling of the short-ranged part 
of certain pair potentials between the ions is performed 
to simulate further covalence effects in the calculation 
in such a way that the energy-minimized structure is as 
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close as possible to the experimental onei^. Structure 
optimization and energy minimization is very important 
for a reliable calculation of the phonon dynamics through 
the dynamical matrix. 

The RIM with the corrections just mentioned then 
serves as an unbiased reference system for the descrip- 
tion of the HTSC's and can be considered as a first ap- 
proximation for the insulating state of these compounds. 
Starting with such an unprejudiced rigid reference sys- 
tem non-rigid electronic polarization processes are intro- 
duced in form of more or less localized electronic charge- 
fluctuations (CF's) at the outer shells of the ions. Es- 
pecially in the metallic state of the HTSC's the latter 
dominate the nonlocal contribution of the electronic den- 
sity response and the EPI and are particularly impor- 
tant in the CuO planes. In addition, anisotropic dipolc- 
fluctuations (DF's) are admitted in our approacb&iii, 
which prove to be specificly of interest for the ions in 
the ionic layers mediating the dielectric coupling and for 
the polar modes. Thus, the basic variable of our model 
is the ionic density which is given in the perturbed state 
by 



p a (r,Qx, Pa ) = p° a (r)+J2QxP\ F (r)+Pa-rp%(r). (1) 



ionic case. Finally, applying the pair-potential approxi- 
mation we get for the total energy: 



a,a (a.a)^(b,/3) 

(2) 

The energy E depends on both the configuration of the 
ions {R} and the electronic (charge) degrees of freedom 
(EDF) {C} of the charge density, i.e. {Q\} and {p a } in 
Eq. (JTJ. E* are the energies of the single ions, a, b 
denote the elementary cells and a, (3 the corresponding 
sublattices. The second term in Eq. @ is the interaction 
energy of the system, expressed in terms of anisotropic 
pair-interactions $ Q /3. Both E* and $ Q/ 3 in general de- 
pend upon £ via p a in Eq. QJ. 

The pair potentials in Eq. (J^J can be seperated into 
long-ranged Coulomb contributions and short-ranged 
terms as follows: 
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p° a is the density of the unperturbed ion, as used in 
the RIM, localized at the sublattice a of the crystal and 
moving rigidly with the latter under displacement. The 
Qx and describe the amplitudes and the form-factors 
of the CF's and the last term in Eq. JIJ represents the 
dipolar deformation of an ion a with amplitude (dipolc 
moment) p a and a radial density distribution r de- 
notes the unit vector in the direction of r. The p^ F are 
approximated by a spherical average of the orbital den- 
sities of the ionic shells calculated in LDA taking self- 
interaction effects (SIC) into account. The dipole density 
is obtained from a modified Sternheimer method in 
the framework of LDA-SIC 10 . All SIC-calculations are 
performed for the average spherical shell in the orbital- 
averaged form according to Perdew and Zunger— . For 
the correlation part of the energy per electron s the 
parametrization given in Ref. ^| has been used. 

The total energy of the crystal is obtained by assum- 
ing that the density can be approximated by a superpo- 
sition of overlapping densities p a . The p° a in Eq. (JIJ 
are also calculated within LDA-SIC taking environment 
effects, via a Watson sphere potential and the calculated 
static effective charges of the ions into account. The 
Watson sphere method is only used for the oxygen ions 
and the depth of the Watson sphere potential is set as 
the Madelung potential at the corresponding site. Such 
an approximation holds well in the HTSC'si2ii4. As a 
general rule, partial covalence reduces the amplitude of 
the static effective charges in mixed ionic-covalent com- 
pounds like the HTSC's, because the charge transfer from 
the cations to the anions is not complete as in the entirely 



$ Q/3 (R, C) =K a Up(R, + K p U a (R, C) 

+ iy a/3 (R,C) + G a/3 (R,C). (4) 

The first term in Eq. @ describes the long-ranged ion- 
ion, the second the dipole-ion and the third and fourth 
term the dipole-dipole interaction. Z a and Zp are the 
variable charges of the ions in case the CF's are excited. 
The latter reduce to the ionic charges for rigid ions. K a 
and Kfj are the charges of the ion cores. The remain- 
ing term in Eq. @ given in represents the short- 
ranged interactions. These short-ranged contributions to 
the pair potentials are expressed by the following inte- 
grals: 



U a (R,Q 
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(5) 



W aP (R,0 = d 3 r dV[p Q (r,C)Mr',C)x 



1 1 (r + rQ-R V 

[ \r-r'-K\ -R- & fi> (6) 

G a/3 (R,C)= J d 3 r[p a p(r,()e(p a p(r,0) 

-p a (v,Oe( Pa (r,()) (7) 
- P/3 (r-R,C)e(Mr-R,C))], 



with 



Pa/3(r,C) =Pa(r,C)+P/3(r-R,C)- 



(8) 
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K a Up(R, C) yields the short-ranged contribution of the 
interaction between the core a and the density pp accord- 
ing to Eq. 0. W a p(R.,Q represents the short-ranged 
Coulomb contribution of the interaction of the density 
p a with the density pp. G a p(R.,Q is the sum of the 
kinetic one-particle- and the exchange-correlation (XC) 
contribution of the interaction between the two ionslfi. 
The short-ranged part of the potentials and the various 
coupling coefficients are calculated numerically for a set 
of disctances R between the ions. The corresponding re- 
sults are than described by an analytical function of the 
form: 



f(R) = ±exp 



(9) 



a, j3 and 7 in Eq. ^ are fit parameters. 
From the adiabatic condition 



dE(R,Q 



= 



(10) 



an expression for the atomic force constants, and accord- 
ingly the dynamical matrix in harmonic approximation 
can be derived: 



t<f(q)= C(q) 



RIM 



j2[Bn<i)r[c- i (c l )] i 



s- ,/3 (q)- 



(11) 

The first term on the right hand side denotes the con- 
tribution from the RIM. M a , Mp are the masses of the 
ions and q is a wave vector from the first Brillouin zone 
(BZ). 

The quantities B(q) and C(q) in Eq. (|ll|l represent 
the Fourier transforms of the electronic coupling coeffi- 
cients as calculated from the energy in Eq. (J2J, or the 
pair potentials in Eqs. @-||HJl, respectivly. 
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k denotes the EDF (CF and DF in the present model, 
see Eq. Q) in an elementary cell. The B coefficients 
describe the coupling between the EDF and the dis- 
placed ions (bare electron-phonon coupling), and the co- 
efficients C determine the interaction between the EDF. 
The phonon frequencies w ff (q) and the corresponding 
eigenvectors e Q (qcr) of the modes (qer) are obtained from 
the secular equation for the dynamical matrix in Eq. 
dJ, i.e. 



^(q)e^(q) 



^ 2 (q)e?(q). 



(14) 



The Eqs. (|11|I - H14|) are generally valid and, in particu- 
lar, are independent of the specific model for the decom- 
position of the perturbed density in Eq. and the pair 
approximation Eq. (j2J) for the energy. The lenghty de- 
tails of the calculation of the coupling coefficients B and 
C can not be reviewed in this paper. They are presented 
111 Ref. [T3 In this context we remark that the coupling 
matrix C KK '(q) of the EDF-EDF interaction, whose in- 
verse appears in Eq. I|ll|) for the dynamical matrix, can 
be written in matrix notation as 



C = IT 



V. 



(15) 



II 1 contains the kinetic part to the interaction C 
and V the Hartree and exchange-correlation contribu- 
tion. C _1 needed for the dynamical matrix and the EPI 
is closely related to the (linear) density response function 
(matrix) and to the inverse dielectric function (matrix) 
e , respectively. 

Only very few attempts have been made to calculate 
the phonon dispersion and the EPI of the HTSC's using 
the linear response method in form of density functional 
perturbation theory (DFPT) within LDA 15 i 16 i 17 . These 
calculations correspond to calculating II and V in DFT- 
LDA and for the metallic state only. On the other hand, 
in our microscopic modeling DFT-LDA-SIC calculations 
are performed for the various densities in Eq. Q in order 
to obtain the coupling coefficients B and V. Including 
SIC is particularly important for localized orbitals like 
Cu3d in the HTSC's. SIC as a correction for a single 
particle term is not a correlation effect, which per defini- 
tion cannot be described in a single particle theory, but 
SIC is important for contracting in particular the local- 
ized Cu3d orbitals. The latter are correlated orbitals and 
so we use sometimes in this paper the phrase " correlated" 
Cu3d states or orbitals, respectively. Our theoretical re- 
sults for the phonon dispersion, which compare well with 
the experiments^, demonstrate that the approximative 
calculation of the coupling coefficients in our approach 
is sufficient, even for the localized Cu3d states, see also 
section IIIII Written in matrix notation we get for the 
density response matrix the relation 



u(i + vuy 



e = l + VTI. (16) 



The CF-CF submatrix of the matrix II can approxi- 
matively be calculated for the metallic (but not for the 
undoped and underdoped) state of the HTSC's from a 
TBA of a single particle electronic bandstructure. In 
this case the electronic polarizability n in tight-binding 
representation reads: 



n KK ,(q,^ = o) = -- ]T 



/„'(k + q)-/„(k) 



x [C* n (k)C Kn ,(k + q)] [C:, n (k)C K/n/ (k + q)]* 



(17) 
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/, E and C in Eq. (|17fl are the occupation numbers, the 
single-particle energies and the expansion coefhcientes of 
the Bloch-functions in terms of tight-binding functions. 

The self-consistent change of an EDF on an ion in- 
duced by a phonon mode (qcr) with frequency kv(q) and 
eigenvector e Q (qcr) can be derived in the form 



-^X KQ (q)u Q (qa) 



iqR a 



with the displacement of the ions 



(18) 



u*(qa) 



2M a ^(q) 



1/2 



e Q (qcr)e lqRa = u Q (qcr)e lqR1 \ 
(19) 

The self-consistent response per unit displacement of 
the EDF in Eq. (|18fl is calculated in linear response the- 
ory as: 



X(q) = n(q) £ - 1 (q)B(q) = C- 1 (q)B(q). 



(20) 



A measure of the strength of the EPI for a certain 
phonon mode (qcr) is provided by the change of the self- 
consistent potential in the crystal felt by an electron at 
same space point r in this mode, i.e. 6V e s(r,qa). Aver- 
aging this quantity with the corresponding density form 
factor p K (r — R*) at the EDF located at R*, we obtain 



6V*(qa)= J dVp K (v-R*)6V cS (r,qa). 



(21) 



This gives a measure for the strength of the EPI in the 
mode (qcr) mediated by the EDF considered. For an ex- 
pression of 6V* (qcr) in terms of the coupling coefficients 
in Eqs. (|12|) . (|13|) . see Ref. 01 From our calculations for 
LaCuO large values for <5V,f (qcr) are found, in particular 
for the phonon anomalies and even larger for the non- 
adiabatic c-axis phonons in the metallic phase, mixing 
with the plasmon. 

The generalization for the quantity II in Eqs. (|15fl . 
(jl6|l needed for the kinetic part of the charge response 
in the non-adiabatic regime, where dynamical screening 
effects must be considered, can be achieved by adding 
[fvjj + irj) to the differences of the single-particle ener- 
gies in the denominator of the expression for II in Eq. 
(|17J) . Other possible non-adiabatic contributions to C re- 
lated to dynamical exchange-correlation effects and the 
phonons themselves are beyond the scope of the present 
model. The coupled-mode frequencies of the phonons 
and the plasmons must be determined self-consistently 
from the secular equation (|llfl for the dynamical matrix 
which now contains the frequency to implicitly via II in 
the response function C _1 . Analogously, the dependence 
on the frequency is transferred to the quantity X in Eq. 
l|2U)l and thus to 8( K and 5V K in Eqs. and (|2T)) . 



respectively. Such a non-adiabatic approach is necessary 
for a description of the interlayer phonons and the charge- 
response within a small region around the c-axis^. 



III. RESULTS AND DISCUSSION 

A. A modeling of the electronic charge response 
and corresponding phonon calculations 

An important question remains as to how to discrim- 
inate between the charge response of the metallic- and 
the insulating state of the HTSC's. The latter cannot 
be obtained within the LDA and a realistic quantita- 
tive description to calculate the irreducible polarization 
part II KK '(q) for the HTSC's is not available. However, 
a general criterion follows from the different analytic be- 
haviour of the polarizability in the long-wavelength limit 
(q — > 0) in both phases^. In the metallic phase the elec- 
tronic partial density of states (PDOS) at the Fermi level 
Zk,(e-f) is related to the polarizability matrix for (q — > 0) 
according to 

^n KK ,(q^0) = Z K (e F ), (22) 
and the total density of states at energy e is given by 



Z{e)=Y,Z K {e). 



(23) 



On the other hand, for the insulating state we obtain 
the sum rules 



]Tn KK ,(q^ 0) = O(q) 



(24) 



and 



0{q 2 ). 



(25) 



The sum YIkk' n KK '(q — > 0) is equal to p 2 K with p 
the average density and K the compressibility of the 
electronic system. The latter provides a measure of the 
gap in the electronic spectrum because K vanishes as 
a function of the chemical potential in the gap region. 
Equations J221 , lO an d (23) > (2H) , respectively, can be 
considered as an orbital resolved closed form to describe 
the metal-insulator transition in terms of the polariza- 
bility or the compressibility of the electronic system, re- 
spectively. This is a primary tool to characterize the 
corresponding ground state. Such orbital selective sum 
rules are particularly useful in the proximity to a Mott 
insulating phase where different from the band insulator 
the internal degrees of freedom, orbital (and spin), ap- 
proximatively survive. In order to fulfil the sum rules 
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FIG. 1: Calculated phonon dispersion of LaCuO in the tetragonal structure in the main symmetry directions A ~ (1,0,0), 
E ~ (1, 1, 0) and A ~ (0, 0, 1) in the model for the optimally doped state. The various symbols representing the experimental 
results indicate different irreducible representations (ID's). Open symbols are for Lai.gSro.iCuO^— and solid symbols for 
Lai.ssSro.isCuOA The arrangement of the panels from left to right according to the different ID's is as follows: 
| Ai | A 2 (- ■ o), A 4 (- □) | A 3 | Ei | E 2 (- ■ ., o), E 4 (- □) | E 3 | Ai(-, □), A 2 (- ■ ■, o) | A 3 



above, in case of the insulating state, in contrast to the 
metallic state, off-diagonal elements of the polarizability 
matrix describing nonlocal polarization processes neces- 
sarily must occur and interfere in order to correlate the 
charge response in such a way that Eqs. (|24p. i|25|) are 
satisfied. 

Physically this is related to the fact that in an insu- 
lator a perturbation, e. g. a change of the electron-ion 
potential, is only incompletely screened. Consequently, 
the selfconsistent change of the potential at the orbital 
degree of freedom n' nonlocally contributes to the CF's 
of the orbital k ^ n' in the unit cell. This is quite differ- 
ent from the metallic case where the diagonal elements 
of dominate and off-diagonal elements can be ne- 
glected in most cases. As mentioned in the Introduction, 
the model for the polarizability of the metallic state is 
based on LDA calculations which, however, have to be 
modified to correct for the overestimation of the c-axis 
coupling in LDA-based investigations 3 . So, a model tak- 
ing into account Cu3d, 4s, 4p, and O xy 2p orbitals in the 
CuO plane results 3 . In contrast to the Cu4s orbital the 
effect of the Cu4p orbital on the phonon dynamics is 
only weak, so we have neglected the latter in the present 
calculations and enhanced the contribution of Tl{CuAs) 
slightly. Besides the CF's also anisotropic DF's are taken 
into account. 

The resulting calculated phonon dispersion is shown in 
Fig. ^ We get an overall agreement with the experimen- 
tal dispersion curves. A few branches show larger devia- 
tions between experiment and theory. However, in par- 



ticular the important phonon anomalies (high-frequency 
oxygen bond-stretching modes), which are a major theme 
in the recent literature and the present paper, too, are 
well described by our theory, see FigQ] below. This also 
holds true for the steep Ai branch which results from 
the large anisotropy of the real material and cannot be 
described with pure LDA based calculations leading to 
an overestimation of the coupling along the c-axis, see 
Ref. Freezing-in the unstable E3 tilt mode at the 
X point correctly indicates the experimentally observed 
structural phase transition from the high-temperature 
tetragonal (HTT) to the low-temperature orthorhom- 
bic (LTO) structure. The actual transition, of course, 
cannot be studied in the harmonic approximation be- 
cause the anharmonic contribution to the energy must be 
considered 1 ®. However, from our calculations in Ref. Q 
we find that the soft tilt mode found in HTT is stabilized 
in LTO in harmonic approximation and thus the struc- 
tural transition is essentially driven by the long-ranged 
ionic forces in the material. 

The non-empirical potential-induced breathing (PIB) 
modeliS like our model reflects the strong tendencies to 
ionic forces in the HTSC's. This model also correctly 
predicts, that the tetragonal structure is unstable to the 
A-point tilt mode. However, it can not be applied to 
the metallic phase. Moreover, other unstable phonon 
branches are found in this model. The width of the 
phonon spectrum is considerably overestimated. For ex- 
ample, the PIB result for the planar oxygen breathing 
mode at the A-point in the insulating phase is about 
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FIG. 2: Same as Fig. 0with the calculated results from the model for the underdoped state. 



34 THz in comparsion to the experiment (~ 21.8 THz) 
and to our calculated result of 21.6 THz, see Fig. 01 
The PIB results are comparable to our calculations per- 
formed with the ab-initio rigid ion model using nominal 
ionic charges^, neglecting the covalence correction for the 
ionic charges and the pair potential. Taking the latter ef- 
fects into account, as mentioned in Section^ the correct 
width of the spectrum is obtained and, with the excep- 
tion of the tilt mode at X, all unstable phonon branches 
disappear. So, we obtain a suitable reference system 
for the HTSC which allows for the investigation of the 
additional nonlocal, non-rigid screening effects in terms 
of CF's and DF's. 

Empirical shell model calculations for the lattice dy- 
namics of several cuprates have been performed in Ref. 
I20L The fitting to the experimental results is quite satis- 
factorily emphasizing again the importance of the ionic 
forces. However, in the fitting the harmonic X-point tilt- 
mode was assumed to be stable, which is at least ques- 
tionable, because of the low temperature orthorhombic 
structure. For a description of the metallic phase free- 
carrier screening of the homogeneous electron gas is addi- 
tional assumed, but the softening of the phonon anoma- 
lies can not be accounted for by this type of screening. 
On the other hand, our calculations, see Fig. lead to 
very good results for the phonon anomalies. Accordingly, 
we can relate the anomalous softening to the characteris- 
tic ionic, inhomogenous electronic structure of the HTSC, 
leading to strong nonlocal EPI in terms of localized CF's 
in the outer shells of the ions controlled by the coupling 
coefficient C in Eq. (|TB|l 

In the model for the insulating phase in order to simu- 
late stronger correlation effects as in the metal Cu4s or- 
bitals are neglected and only the CF's of the Cu3d or- 



bitals at the Cu and 02p orbitals at the O xy are con- 
sidered as on site CF's and the polarizability is adjusted 
so as to obey the sum rules from Eq. I|24|) . The cal- 
culated results for the phonon dispersion can be found 
in Ref. 0- Similarity as for the metallic state a good 
agreement with the experiments is obtained. In par- 
ticular, a significant reduction of the phonon anomalies 
as compared to the optimally doped material is well re- 
produced and attributed to the incompressibility of the 
electronic state. Including the screening contribution via 
anisotropic DF's, the longitudinal and transverse optical 
modes of Am and E u symmetry are well described. 

An ab initio calculation of the polarizability n for the 
HTSC's, which besides strong electron correlations also 
should include the strong nonlocal EPI found in our cal- 
culations in the self-energy of the electrons, seems not 
possible in the near future. In particular this holds for the 
insulating and underdoped state. However, modeling n 
according to the the sum rules l|22|) - (|25() guarantees that 
the constraints on the charge response in the different 
phases are fulfilled and, moreover, allows for the possibil- 
ity to theoretically design condensed matter states with 
properties between that of a conventional metal and an 
insulator, respectively. This seems necessary to describe 
the physics in the underdoped cuprates where a strange 
metallic state does appear. Our modeling of the cor- 
relation phenomenon and the delocalization-localization 
transition from the metallic state of the HTSC's via the 
underdoped state to the undoped (insulating) state has 
been achieved with help of the sum rules by consecutive 
orbital selective compressibility-incompressibility transi- 
tions, see Ref. 0, and for a detailed review Ref. Q]. Start- 
ing from the metallic side, representing in our model the 
doping range from optimally doped to overdoped, all the 
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admitted orbitals, i.e. Cu3d, 4s and 02p, are taken to be 
compressible, metallic consistent with an emerging quasi- 
particle description and a well defined Fermi surface. The 
corresponding sum rules appropriate for such a situation 
are given by Eqs. JUJ), with all the Z k (ef) differ- 
ent from zero. Our calculation of the orbital occupation 
numbers^ from a TBA of the first principles electronic 
(metallic) bandstructure of LaCuO 21 indicates that be- 
sides the Cu3d in particular the more delocalized Cu4s 
state must be included as an additional EDF at the Cu 
site. This, as already mentioned, is different from our 
modeling of the insulating - and underdoped state where 
we have neglected a contribution from the Cu4s orbitals 
and only used the localized Cu3d - and the 02p orbitals 
as EDF's in order to simulate stronger correlation effects 
in these states. 

Approaching the delocalization-localization transition 
from the metallic region when p-doping is decreased first 
the Cu3d component of the wavcfunction is admitted in 
our model to become incompressible, insulator-like in the 
underdoped regime by assuming that the corresponding 
single particle density of states Z K (e) is suppressed at the 
Fermi level for the correlated localized Cu3d orbitals, be- 
cause of the cost in energy from hopping of the charge 
carrier to the Cu sites, but not so for the more delo- 
calized 02p orbitals at the O xy -sublattices in the CuO 
plane where the holes are predominantly injected to in 
p-type cuprates. So, we have an orbital selective com- 
pressible, metallic charge response for the 02p orbitals 
with a renormalizcd PDOS Z K (e-p) and the holes tend 
to accumulate in the region of suppressed antiferromag- 
netism. Altogether, we have a loss in the density of states 
(pseudogap, correlation gap) at the Fermi level. Thus, 
the quasiparticle picture (Z k (ef ^ 0) for all n) consis- 
tent with our modeling of the metallic state for optimally 
to overdoping no longer holds in the underdoped state. 
The related compartmentalization of phase space can be 
expected to be the reason for the strange metallic be- 
haviour observed in the underdoped cuprates with an 
ill defined Fermi surface. The corresponding electronic 
state for this model (pseudogap model, PGM) express- 
ing the dichotomy between a localized and delocalized 
component has been defined by the following sum rule^, 



charge response of the Cu3d orbital at the Cu sublattice 
and a compressible, metallic charge response of the 02p 
orbitals at the O xy -sublattice in the CuO plane is dis- 
played in Fig. Inland compared with experimental inelas- 
tic neutron results for underdoped LaCuO. The agree- 
ment of the calculations, which also include anisotropic 
DF's besides the CF's of Cu3d- and 02p type, is good 
and of similar quality as the calculated results for opti- 
mally doped state in Fig. ^ an d the insulating state in 
Ref. 0- Comparing the calculations in Figs. ^ 121 essen- 
tial differences are only found for the phonon anomalies 
(highest Ai and Si-branch, respectively) which are less 
pronounced in the underdoped state, because the partial 
" correlation gap" introduced by both, the incompressible 
Cu3d orbitals and the missing of the compressible Cu4s 
orbitals leads according to Eq. (|18|l to reduced insulator- 
like CF's at the Cu ion. Compared with the insulating 
state, however, the anomalies are increased in the under- 
doped phase in particular for the A i -branch. A more de- 
tailed discussion of the phonon anomalies is given below. 
Another difference of the theoretical results as compared 
with those of the insulating state given in Ref. □ is of 
qualitative nature and is related to the axial polarized 
Ai-branches. In the modeling of the underdoped state 
we have partially incompressible regions, however, the 
total compressibility is never zero in this state, because 
a real space organization of the low lying charge excita- 
tions is allowed and a metallic (gapless) charge transport 
is enabled via the hole-doped 02p orbitals at the oxy- 
gen network in the CuO plane. As a consequence the 
A^ u mode splittings are closed due to metallic screening, 
in particular the very large "ferroelectric" mode split- 
ting, observed and calculated in the insulating stated 
and the Ai-branch with the steep dispersion characteris- 
tic for metallic behaviour, being absent in the insulating 
state, reappears as is visible from Fig. [21 Altogether, 
we conclude that our modeling of the insulating, under- 
doped and optimally doped state of LaCuO is consistent 
with the measured phonon dispersion. 



B. Qualitative discussion of the modeling of the 
electronic state in the cuprates 



' \Z K {s F ) 



Cu3d 
02p. 



(26) 



We realize from this equation, that nonlocal polariza- 
tion processes are necessary to depress the density of 
states at £f via the Cu3d component and correspondingly 
the low energy single particle excitations. Finally, near 
half filling in addition to Cu3d, the 02p component of 
the electronic state has been modeled as incompressible 
too, and the charge-transfer insulator is formed accord- 
ing to the sum rules (|24ll . (|25|l for both the Cu3d and 
02p orbitals. 

The phonon dispersion curves for our model of the 
underdoped state with an incompressible insulator-like, 



We now discuss some qualitative aspects of our mo- 
deling of the electronic state of the cuprates suggested 
by the phonon calculations. In the model for the under- 
doped state a metallic mobility behaviour is mediated 
via the 02p component of the wavefunction localized at 
the oxygen network of the CuO plane where the holes 
are injected too, but blocked along the Cu-0 links be- 
cause of the large on-site repulsion of the Cu3d compo- 
nent. Thus, by this partial ordering we have hole rich and 
hole poor regions of the conducting carriers on a micro- 
scopic scale in a translational invariant way ultimately 
driven by competing energy contributions. Note, that 
the various inhomogeneous selforganized charge and spin 
ordered configurations discussed in the litertatur o 22 i 23 i 24 
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formally can be considered as certain " defect structures" 
(e.g. stripes) of our model if the spin-degrees of free- 
dom are explicitly considered. In case such structures 
are experimentally observed they express the tendency 
to restore an insulating Mott state in certain extended 
(incompressible) regions of the material. 

From an energetic point of view the ground state 
should simultaneously minimize the zero point kinetic 
energy of the doped charge carriers and the exchange 
energy where both compete which each other. Because 
of the incompressible localized Cu3d states the holes are 
expelled from the antiferromagnetic regions by Coulomb 
correlation. Thus, an independent particle picture where 
an electron feels a time-averaged local density of the other 
electrons with a corresponding local potential (like e.g. in 
LDA) no longer holds. In this way the exchange energy 
is optimized (no spins need to be switched because of the 
restricted motion of the holes) while on the other hand, 
there is a gain in kinetic energy in the hole rich regions by 
derealization of the electronic structure (hopping of the 
holes) via the compressible, metallic 02p component of 
the wavefunction. At the same time the exchange energy 
can be neglected in this region of configuration space if 
the derealization effect is strong enough. Spin degrees 
of freedom and their correlations are not explicitly in- 
cluded in our modeling but implicitly support the latter 
because residual antiferromagnetic correlations can lead 
to an additional lowering of the kinetic energy by vir- 
tually hopping to the Cu3d states. This correlates the 
charge transport in the regions of excess charge to the 
remaining antiferromagnetic spin order in the hole poor 
regions at the Cu sublattice. So, we have a picture that 
can be regarded as resulting from a kind of separation 
of spin and charge, being mutually relevant in different 
regions of the CuO plane. Low energy charge excitations 
occur in the hole rich regions where spin effects are de- 
pressed but are blocked in the hole poor regions, where 
the spin degrees of freedom are most important at low 
energy and high-energy (virtual) CF's induce an antifer- 
romagnetic coupling on the Cu sublattice. 

Then, in our modeling of the metallic state at higher 
doping the Cu3d and the additionally allowed delocalized 
Cu4s component become compressible, metallic. The in- 
compressibility of Cu3d can be thought to be destroyed 
by an increased population of the Cu4s orbital. In Ref. 
l25l it has been shown that the Cu4s admixture to the elec- 
tronic states of the HTSC is relevant and due to a gain 
in correlation energy. The dominant charge transfer re- 
sulting from the atomic correlations is from the localized 
Cu3d to the more delocalized Cu4s orbitals. The calcu- 
lations in Ref. |25| are based on an ab-initio method (local 
ansatz method) that adds correlations as corrections to a 
single-particle self-consistent ground state obtained from 
a Hartree-Fock calculation. An increased hybridization 
of Cu4s with the pda orbitals can be expected to further 
relaxing the kinetic energy when doping is increased. As 
a consequence, the blocking of the metallic charge re- 
sponse at the Cu sites in the underdoped state is lifted 



and metallic behaviour is now also possible via the CuO 
links. Altogether, this seems to be consistent with an 
emerging quasiparticle picture, a progressively weaken- 
ing of the antiferromagnetic fluctuations and the devel- 
opment of a well defined Fermi surface as described by 
the sum rules from Eqs. I|22|) . I|23|l in our approach. The 
growing importance of the Cu4s component in the wave- 
function leading to enhanced delocalized CF's at the Cu 
sites in the screening process is also mapped to a cor- 
responding renormalization of the phonon anomalies as 
compared to the situation where the latter are neglected, 
see Ref. l26land the discussion below. So, from our calcu- 
lation additionally to Cu3d and 02p an increasing contri- 
bution of the delocalized Cu4s component to the metal- 
lic state of the cuprates is important to understand the 
charge response and the corresponding phonon softening 
in the real material at least at higher doping levels. 

The importance of the Cu4s orbital for a realistic de- 
scription of the electronic structure of the HTSC's has 
also been pointed out in Ref. |2?] and Ref. |2^. In the 
latter work, the hopping range has been identified as 
an essential material-dependent parameter and the in- 
tralayer hopping beyond nearest neighbours as well as 
the interlayer hopping proceeds via the Cu4s orbital. It 
ist further concluded that materials with higher Tc ;max 
have larger hopping ranges and in materials with high- 
est Tc jlnax the axial orbital, essentially a hybrid between 
Cu4s, Cu3d 3z 2_ 1 and apical oxygen 2p z , is almost pure 
Cu4s. Moreover, the importance of Cu4s for an accu- 
rate description of the partial charge distributions in the 
HTSC's is pointed out applying the "local ansatz" ab ini- 
tio method, which allows to add correlations as correc- 
tions to a Hartree-Fock ground stated or density func- 
tional theorjiSi, respectively. 

Now we address the question concerning the possible 
nature of a superconducting state of the HTSC's that 
would be consistent with our modeling of the normal 
state. The suppression of the local CF's in our model for 
the underdoped state related to the incompressible Cu3d 
orbitals due to the large on-site repulsion competes with 
the superconductivity because of the enhanced phase 
fluctuations according to the number-phase uncertainty 
relationship ANA(f> > 1. In such a phase disordered 
superconductor pairing without long-range phase coher- 
ence is possible, i.e. there is an intermediate temperature 
range (T > Tc) where electron pairs exist which have not 
condensed. Such a precursor pair formation scenario is 
for example discussed in context with a kinetic energy 
driven pairing mechanism related to certain self orga- 
nized local charge inhomogenities in form of stripes22i2£ 
where pairing originates from strong repulsive interac- 
tions. Pair formation is thought to be realized in form 
of spin singlet pairs resulting in a spin gap by the "spin 
gap proximity" effect because it seems favourable that, 
under appropriate circumstances, the pairs can spread 
out somewhat into the antiferromagnetic neighbourhood 
of a stripe. However, up to now there is no experiment 
that proves the pairing is either kinetic energy driven or 
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due to stripes. Note in this context, that the stripe sce- 
nario has been questioned by recent experiments which 
find that the observed nanoscale disorder in the HTSC's 
is due to dopant disorder—. Moreover, if the possibil- 
ity of precursor pairs of some kind is theoretically ex- 
plored, the strong nonlocal EPI found in our calculations 
expressed by the corresponding large changes of the self- 
consistent potential an electron feels according to Eq. 
Q21[l. in particular for the generic phonon anomalies 31 
and even larger for the non-adiabatic c-axis phononsi^, 
should be taken into account in the electronic self-energy 
and the polarizability besides the short-ranged electron- 
electron repulsion that leads to spin- and charge corre- 
lations in a many body treatment of a preformed pair 
state and the resulting gaps in the normal as well in the 
superconducting state. For a recent review of the impor- 
tant role of the EPI for the normal state properties and 
the pairing in the HTSC's we also refer to Ref. |33- Al- 
together, in the normal state one can expect two types 
of gaps in our modeling for the underdoped HTSC's the 
" correlation gap" and consistent with that modeling pos- 
sibly a " precursor pairing gap" which are both linked by 
the short-ranged correlations introduced through the lo- 
calized, incompressible Cu3d orbitals. While the former 
gap should correspond to the high-energy pseudogap seen 
in ARPES, the latter might be related to the low-energy 
gap found by the same experimental technique 6 . 

Different from the underdoped state in our modeling 
of the optimally doped metallic state all the orbitals be- 
come compressible, metallic with enhanced CF's partic- 
ularly in the Cu4s orbitals. So, we have a qualitative 
change to an electronic ground state with a well defined 
Fermi surface. In such a state a Fermi surface instabil- 
ity with a simultaneous evolution of pairs of quasi parti- 
cles would be consistent with the normal state and can 
lead to a BCS-type ordering. With regard to predom- 
inantly a phase ordering transition in the underdoped 
state optimal doping with a maximal Tq can be thought 
of as defining a gradual crossover to an overdoped state 
with predominantly a BCS-type pairing transition. Ar- 
guments in favour of such a crossover scenario also have 
been presented in Ref. I22H23I 

According to the modeling of the charge response sup- 
ported by our phonon calculations there are higher en- 
ergy CF's related to the localized Cu3d orbital and low 
energy CF's of the more delocalized 02p - and espe- 
cially the Cu4s orbitals and a critical mixing of both 
via a population of Cu4s and a corresponding change 
of the populations of Cu3d and 02p and the related 
hybridization obviously is required to achieve Tc jlnax 
in the optimally doped material. The localization- 
delocalization transition in the normal state expressed 
as an incompressibility-compressibility transition of the 
Cu3d states corresponds to a phase ordering-pairing tran- 
sition in the superconducting state. 



C. Doping dependent anomalous phonon branches 
- Nonlocal dynamic charge inhomogenities 

Next, we report new calculations of the phonon disper- 
sion for the overdoped state of LaCuO which can be com- 
pared with recent experimental results for the phonon 
anomalies^ 3 -. An important message from these calcula- 
tions is a further growing of the importance of the delo- 
calized Cu4s component in the electronic wave function 
which seems consistent with an improved Fermi liquid 
quasi-particle picture, an increased coupling along the 
c-axis and the dying out of the antiferromagnetic corre- 
lations, and an EPI which becomes more local. 

In earlier calculationsi*2& we already have recognized 
the inclusion of Cu4s into the orbital basis as an im- 
portant "tuning parameter" for the derealization of the 
electronic structure of the HTSC's and the strength of the 
phonon anomalies. A growing occupation of the Cu4s or- 
bital and a corresponding growth of the Cu(4s)-element 
of the polarizability matrix leads to an increased soft- 
ening of the anomalous phonon modes. So, enhancing 
n(Cu4s) as compared with its value in the optimally 
doped state leading to the results shown in Fig. ^ gives 
very good agreement with the values measured for the 
anomalies in the overdoped stated and, moreover, pro- 
vides a prediction for the complete phonon dispersion not 
measured so far, see Fig. [3J 

In Fig. 0] we show a comparison of our calcu- 
lated results for the anomalous phonon branches, Ai 
and Si, in LaCuO with the corresponding experimen- 
tal dispersion 3 ^. The strong doping dependence of the 
experimental phonon frequences and thus of the EPI is 
clearly visible and very well accounted for by our mode- 
ling of the electronic state of the HTSC's discussed above. 
We have also included the phonon dispersion as obtained 
for the RIM representing approximatively the ionic local 
part of the charge response and the EPI in order to em- 
phasize the large contribution due to the nonlocal EPI 
and screening in terms of the CF's. The displacement 
patterns of the phonon anomalies, i e. the planar oxy- 
gen breathing mode at the X point, Og , and the oxygen 
half-breathing mode, Ai/2, are displayed in Fig. [SJ For 
symmetry reasons there are only CF's at the silent Cu 
ion allowed in Og and not on the moving O xy ions while 
in Ai/2 in addition to the CF's at the Cu ion also CF's 
can be excited at the silent oxygens. This leads to an ad- 
ditional source of softening in case of Ai/2 as compared 
to Og and the larger renormalization of the frequencies 
found in our calculations and the experiments 3 ^* 3 ^. As 
already mentioned the enhanced softening in the over- 
doped doped state can be attributed to the growing Cu4s 
component in the electronic state and the related CF's 
at low energy. Moreover, the local character of the EPI 
is increased being dominant in a normal metal and su- 
perconductor, with delocalized electrons. 

A word should be said concerning the importance of a 
correct representation of the on-site Coulomb repulsion, 
Ud, related to the localized Cu3d states and their effect 
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(a) Experimental Results (b)Calculated Results 

FIG. 4: (a) Experimental results for the highest Ai and £i branch of I^-^Sr^CuCU for various doping levels—: »x — 0.00, 
Mx = 0.10, ox = 0.15, Ox — 0.30. The lines are a guide to the eye. (b) Calculated results of the phonon branches shown in 
(a) as obtained within our model approach for the electronic state of the HTSC's: — insulating state, ■ ■ ■ underdoped state, - 
- optimally doped state, — • — overdoped state. For comparison, the results for the ab initio RIM are shown (—•••) in order 
to demonstrate the large influence of the nonlocal EPI effects in form of CF's on the dispersion. 
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FIG. 5: Displacement patterns of the high-frequency oxy- 
gen bond-stretching modes of LaCuO. Left: planar breath- 
ing mode (Og), middle: half-breathing mode (Ai/2), right: 
apex-oxygen breathing mode (Of ). 

on the nonlocal EPI by suppressing low energy CF's. The 
onsite-parameters, like Ud for the Cu3d orbitals, are de- 
fined in our modeling using Eq. (H3Jl and taking the EDF 
re, re' to be localized at the same ion a in the crystal. An 
explicit expression for the calculation of these parameters 
is given by the following equations^. 

+ JdVp^T)p,.(T)m[pHr)], (27) 

with 

w{ P ) = = (28) 

The relevant form factors p K , p K < in Eq. ijTH) are de- 
fined by the pf F from Eq. and p° a is the unperturbed 
density of the ion a calculated in LDA-SIC as noted in 
Section|nJ The calculated value for Ud in units of e 2 /aB is 
1.0050. The corresponding value for the more delocalized 
02p and Cu4s states is 0.6708 and 0.3563, respectively. 
The large value for the Cu3d orbitals tends to supress 
CF's at the Cu site according to Eqs. (|T5|> . (fTHf) and 
()20JI , while the latter arc strongly enhanced by the much 
smaller value for the Cu4s orbital. This again points to 
the important role of a Cu4s admixture in the electronic 
state and for the anomalous phonons in the cuprates. 

A characteristic feature of the phonon anomalies is 
their selective sensitivity to Ud- A decrease of Ud en- 
hances the CF's at the Cu ion considerably and leads 
in parallel to a softening of the phonon anomalies while 
the other phonon branches coupling to CF's are virtu- 
ally not affected; see Fig. [S]for the Ai/2 anomaly, where 
we have decreased Ud from its calculated value which 
by itself leads to the good quantitative agreement with 
the experiments in Fig. [21 A similar calculation for the 
Ei branches even leads to an unstable planar breath- 
ing mode Og being the highest mode in the insulating 
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FIG. 6: Demonstration of the selective sensibility of the 
anomalous highest Ai branch with respect to the on-site 
Coulomb repulsion Ud of the localized Cu3d orbitals. Only 
Ai branches coupling to CF's are shown for the model of un- 
derdoped state of LaCuO. Using the calculated value for Ud 
(DFT-LDA-SIC, left panel) good agreement with the exper- 
imental data (Fig. |5J is obtained. In the remaining panels 
from left to right we arbitrarily have decreased Ud by the fac- 
tor |, i and |, respectively. For ^Ud the frequency of the 
Ai/2 anomaly (~4.8 THz) is decreased towards an instabil- 
ity. In a corresponding calculation for the Ei branches Ob 
even becomes unstable for \Ud- 



and underdoped state (other phonon branches coupling 
to CF's are only weakly modified). So, we conclude that 
our approximative calculation of Ud within DFT-LDA- 
SIC is sufficient to describe the correlation effect related 
to Ud on the phonon dynamics. 

Local charge transfer from Cu3d into Cu4s orbitals due 
to atomic correlations has been shown in Ref. |23 to lead 
to a strong reduction (screening) of Ud for effectiv model 
Hamiltonian calculations. The influence on the phonon 
anomalies in LaCuO of such a screening effect for Ud al- 
ready has been investigated in Ref. 0j. In this work a 
model with a reduced set of EDF is used and CF's are 
only allowed for the Cu3d and the O xy 2p states of the 
ions in the CuO plane. The contribution of more delocal- 
ized states, like Cu4s or Cu4p, is simulated in this model 
by assuming a renormalizing (screening) of the on-site 
interaction Ud- Within such a procedure, the anomalous 
softening for Ai/2 and Og now seen in the experiments 
has been predicted and related to the changes of the effec- 
tive repulsive short-range Coulomb interaction at the Cu 
ion as an important physical parameter of the cuprates. 

From an ab-initio point of view the definite value for 
an effectiv Ud parameter which should be used in mod- 
els with a reduced set of electronic degrees of freedom, 
like the Hubbard models, which is widely used in the 
cuprates is by no means straightforward. Quite gener- 
ally, a screening of the Coulomb interaction can be dis- 
cussed within a renormalization scheme of the electronic 
density response as proposed iniLi&Si, leading for the 
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inverse dielectric matrix e 1 to the following identity 

£- 1 =e- 1 -e-\ (29) 

with 

£~ l = (1 + uA)" 1 , e" 1 = (1 + vfiy 1 , w = e" 1 i;. 

(30) 

This general factorization of the screening follows from 
an additive decomposition of the total polarizability II 
into a suitable chosen renormalizing contribution (II) and 
a remaining part , A, that contains all the polarization 
processes not included in II, i.e. 

II = n + A. (31) 

Applying Eq. (|29() to the bare Coulomb interaction, v, 
we obtain the screened interaction 

V sc =s- l v= (1 + vA)~ 1 v. (32) 

In (|31|l we can for example identify A with that part 
of II which contains only 3d-3d transitions and Tl with 
the remaining part of the total polarizability, i.e. by the 
other (renormalizing) electronic orbitals considered. Us- 
ing this specific realization of the decompostion (|31|l , we 
can extract from Eq. I|32|) that the interaction between 
the 3d electrons is provided by the interaction v defined 
by the screening processes of the other orbitals taking 
into account in II. 

In general the interaction between the 3d electrons, 
v, becomes frequency-dependent (retarded) via the 
frequency-dependence of the polarizability II, which 
could be calculated e.g. at the RPA level, like in our 
treatment of the non-adiabatic regime. Summarizing, 
the interaction in a model reduced to the 3d degrees 
of freedom should be frequency- dependent in contrast to 
the static t/d-parameter assumed intuitively in Hubbard- 
like Hamiltonian approaches. From our considerations 
a static on-site interaction Ud should be identified with 
the partially renormalized interaction v in the low-energy 
limit, i.e. Ud — v(oj — * 0). In this context we refer to 
recent work where it is found that a static value of Ud 
may not be the most appropriate one to use in effective 
model calculations^. 

In summary, from our calculations we find that the 
phonon anomalies and the strong nonlocal EPI can be 
understood in terms of excited CF's in the outer shells of 
the ions, which are controlled by the kinetic contribution 
II -1 to the electronic coupling coefficient C in Eq. 115|) 
together with the contribution V, i.e. the Hartree- and 
exchange-correlation part of C. The interplay between 
strong correlation effects related to the localized Cu3d 
orbitals with their tendency to suppress CF's with low 
energy via V and the contribution of the more extended 



orbitals 02p and Cu4s with a smaller (calculated) on- 
site repulsion which tends to enhance CF's is an essential 
feature for an understanding of the charge response and 
its influence on the phonons in a quantitative manner. 

We have made explicit the charge response by 
calculating the phonon-induced charge density 
redistributioni2iS& for the phonon anomalies. Ac- 
cording to these investigations the nonlocal EPI 
generates dynamic charge ordering via CF's in form of 
localized stripes of alternating sign in the CuO plane 
leading to an effective reduction of the energy visible 
in the anomalous softening. For the Ai/2 mode the 
charge stripes point along the x- or y-axis, respectively, 
and for Og along the diagonals in the CuO plane. 
These dynamical charge inhomogeneities would not be 
present in systems with local EPI, for example in a 
high density homogeneous electron gas prevailing in 
conventional metals. On the other hand, certain charge 
inhomogeneities observed in the HTSC's because of the 
reduced screening will couple nonlocally to the lattice 
and additional softening, compared with the system free 
of inhomogeneities, of the oxygen-bond-stretching- like 
modes (OBSM) can be expected due to the sensitivity 
of the latter to charge doping as seen in Fig. For 
example, in the presence of a periodic charge order at a 
certain wavevector compatible with its period, anoma- 
lous softening and related mode splitting according to 
the broken rotational and translational symmetry should 
occur. Also non-periodic, local charge inhomogeneities 
could lead to anomalous local OBSM via non-local 
EPI. These theoretical expectations are supported by 
recent experiments^. Here strong phonon softening 
of the OBSM has been observed in HTSC's at doping 
levels associated with static stripes at a wave vector 
corresponding to the charge order. From the calculated 
results of the anomalies in Fig. 0] for the system with 
unbroken symmetry it can be expected that the strength 
of the anomalies in the inhomogeneous case sensitively 
should depend on the possibility of the inhomogeneous 
structure (e. g. stripe) to be incompressible, insulating 
or compressible, metallic, respectively. 

It is interesting to note that the electron correlation 
which result in an antiferromagnetic spin correlation in 
the HTSC's should favour the CF's in the OBSM and 
enhance screening against the tendency of a suppression 
of the CF's by the large on-site repulsion at the Cu, be- 
cause if the spins in neighbouring Cu orbitals were to be 
parallel, they would not both to be able to transfer to 
the neighbouring oxygens or between the Cu ions them- 
selves. On the other hand, if they are antiparallel they 
can do. Thus the phonon induced CF's depend implicitly 
on the spin-degrees which may enhance the former and 
the EPI in an antiferromagnetic neighbourhood as com- 
pared to a paramagnetic or a ferromagnetic one. Vice 
versa the nonlocal EPI of CF-type is expected to gener- 
ate antiferromagnetic spin-fluctuations via the phonon- 
induced transferred charge between the Cu ions. This 
discussion around the phonon anomalies provides an ex- 
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ample to illustrate, how coupled lattice-charge- and spin 
degrees of freedom could act synergetically for pairing in 
the cuprates in a doping dependent way. 

Besides Ai/2 and Og essentially polarized in the CuO 
plane another OBSM, i.e. the apex-oxygen (O z ) bond- 
stretching mode at the Z point, (Of), polarized perpen- 
dicular to the CuO plane, (Fig. EJ, is of special inter- 
est for the cuprates due to the long-range Coulomb in- 
teraction in these compounds. Here we have the situ- 
ation (similar as for the Laf mode) that the displace- 
ment of the ions (O z , La) in the ionic layers generates 
changes of the potential felt by the electrons in the CuO 
plane which themselves are responsible for the super- 
conductivity. Such nonlocal coupling effects are an ex- 
pression of the strong component of the ionic binding 
along the c-axis in the HTSC's, i.e. these long-range 
Coulomb coupling effects are very special to this class 
of materials, and would not be possible in a convential 
metal or superconductor because of local-screening by 
a high-density electron gas. From the experiments in 
LaCuO^^ we get the very remarkable result that Of 
is strongly renormalized when going from the insulat- 
ing parent compound («17 THz) to the optimally doped 
material («11.5 THz). However, in the optimally doped 
probe a very large linewidth of about 4 THz appears that 
makes it difficult to localize the exact position of Of in 
the metallic phase. The massive line broadening has been 
interpreted within our non-adiabatic calculations^ lead- 
ing to a phonon-plasmon scenario for modes propagat- 
ing in a very small region parallel to the c-axis. From 
the calculations^* we find a small region nearby the c- 
axis with a non-adiabatic, insulator-like charge response 
crossing over to a coherent adiabatic, metallic response 
outside this region. Moreover, an interconnection be- 
tween the in-plane response and the size of the non- 
adiabatic inter-plane response is pointed out. So the elec- 
tronic properties in the CuO plane take influence on the 
charge response perpendicular to the plane nearby the 
c-axis. 

The strong softening in the metallic state outside the 
small region of non-adiabatic charge response (for details 
see Ref. Ill ll) can physically be understood by compar- 
ing the phonon induced charge rearrangements for Of as 
calculated in the insulating and metallic state (Figs. 
IB}, respectively. For the frequencies we find a renormal- 
ization of about 3 THz from 17.13 THz in the insulator 
to 14.17 THz in the metal, compare with the adiabatic 
phonon dispersion in Fig. 1 where Of is the second high- 
est frequency of the Ai-branches at the Z point. A large 
EPI for the Of mode in the metallic state of LaCuO 
has also been calculated in Ref. |38| and attributed to a 
combination of weak screening along the c-axis and non- 
local Madelung-like interactions. The importance of such 
long-ranged ionic contributions to coupling also has been 
pointed out in Refs. l39ll40ll44Tl45l 

As can be extracted from Fig. [SJfor Of the apex oxy- 
gens move in phase against the CuO layers. So, because 
of the weak screening one can expect this vibration to in- 
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FIG. 7: Contour plot of the nonlocal part of the phonon- 
induced charge density redistribution for the Of mode (Fig. 
I^J in the model for the insulating state of LaCuO. The units 
are 10 _4 e 2 /a|. Full lines ( — ) mean that electrons are ac- 
cumulated in the corresponding region of space and the lines 
(—•■•) represent regions where the electrons are pushed away, 
resulting in an intra-layer charge transfer. 

duce CF's in the CuO planes. These CF's, however, are 
qualitatively different in the insulating and the metallic 
state, respectively, because of the gap for charge exci- 
tations in the insulator. This helps to understand the 
anomalous softening of the mode during the insulator- 
metal transition. As has been shown in Ref. |5| the CF's 
<5£ K (qer) from Eq. I|18fl are constrained in the insulator 
with CF's allowed at the Cu and O xy sublattices accord- 
ing to the following sum rule 

^<5C K (Aa) = 0. (33) 

K 

The sum over k denotes the CF's in the CuO layer. 
A ~ (0, 0, 1), so Eq. (|3l5j) particularly holds at the Z point 
of the Brillouin-zone. In contrast to the constraint ex- 
pressed by Eq. (|33|l for the insulating state, which means 
that local charge neutrality of the cell is maintained un- 
der a perturbation due to Of, no such a restriction is 
present for the metallic state. Unlike for LaCuO in the 
insulating state, where only intra-layer charge rearrange- 
ments according to the sum rule Eq. I)33|) are allowed 
(Fig. 0|, in the metallic state Of induces CF's at Cu 
and O xy that have the same sign in the whole CuO layer 
(Fig. |SJ|. This finally leads to CF's of alternating sign 
in consecutive CuO layers, i.e. to an inter-plane charge 
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FIG. 8: Same as Fig. [7|but calculated with the model repre- 
senting the optimally doped metallic state, resulting in CF's 
of the same sign in a CuO layer and a corresponding inter- 
layer charge transfer. 



transfer in the adiabatic regime that provides an effective 
screening mechanism for the long-range Coulomb inter- 
actions and generates the anomalous mode softening of 
Of during the insulator-metal transition. 

D. Phonon-density of states 

In the last topic of this paper we present calculations 
for the phonon-density of states of LaCuO according to 
our modeling of the insulating-, under-, optimally- and 
overdoped state of the HTSC's. Experimental studies by 
inelastic neutron spectroscopy of the phonon-density of 
states of LaCuO report unusual shifts for high-frequency 
OBSM during the insulator-metal transition indicating 
a strong EPI 41 . Similar results are found for LaCuO in 
Ref. l42l where besides a shift of the high-frequency part 
of the spectrum to lower frequencies the development 
of " new" " oxygen lattice vibrations centered in the opti- 
mally doped probe around 70 meV are observed near the 
doping induced insulator-metal transition and attributed 
to the mixing of these vibrations with metallic CF's. 

Our calculations for the phonon-density of states dis- 
played in Figs. EH EU allow for a detailed study of the 
changes of the spectrum during the insulator-metal tran- 
sition. With the help of the calculated atom-resolved 
phonon-density of states in Fig. I1UI and the phonon dis- 
persion the characteristic modes indicated in Fig. [!|]can 
be identified. The peak denoted as A is caused by vibra- 



tions of the apex-oxygen, O z . B is related to the high- 
frequency oxygen vibrations, E u , at the T point polar- 
ized in the CuO plane. C is the planar oxygen breathing 
mode, Og , and D results from the oxygen half-breathing 
mode (Ai/2 mode). In agreement with the experiments 
we find globally a shift of the high-frequency part of 
the phonon-density of states and also the development 
of spectral weight (peak D, see also Fig. EU( C )) passing 
from the insulating to the metallic state which we can 
assign to the Ai/2 anomaly at 69 meV in our calcula- 
tions for the optimally doped phase. Thus, the "new" 
oxgen lattice vibrations found in the experiments consist 
dominantly of the half-breathing modes, as supposed in 
Ref. ^3 An important change of the spectrum related 
to the insulator-metal transition can be extracted from 
Figs. El El because B(E U ) and C(Og ) are interchanged. 
So we have not just a global softening but also a very 
characteristic rearrangement of the high-frequency part 
of the spectrum, linked to the development of the phonon 
anomalies and the corresponding softening upon doping. 



IV. SUMMARY 

Within our microscopic modeling of the insulating-, 
underdoped-, optimally- and overdoped electronic state 
of the HTSC's we have performed calculations for LaCuO 
of the complete phonon dispersion curves, dynamic 
charge inhomogeneities induced by nonlocal EPI and the 
total-and atom-resolved phonon density of states. Our 
calculations emphasize the role of the doping dependence 
of the generic phonon anomalies and their selective sen- 
sitivity to both localization (correlation) effects provided 
by the Cu3d orbitals and derealization effects via the 
Cu4s states. Thus, these anomalies are suitable probes 
to mirror the localization-delocalization transition of the 
electronic structure of the cuprates in terms of phonons. 
We have found good agreement with the correspond- 
ing measured phonon dispersion and phonon-density of 
states and have presented an interpretation of the re- 
sults. In order to reach these results, both, long-ranged 
Coulomb interactions of ionic origin as well as short- 
ranged repulsiv interactions together with a sufficient set 
of orbital degrees of freedom are neccesary. A purely 
ionic model leads to a considerable overestimation of the 
width of the phonon spectrum and to unstable branches 
besides the tilt mode. Including additionally covalent 
corrections a suitable starting model for the HTSC can 
be found that cures the disadvantages of the pure ionic 
description. Using such a model as an unprejudiced rigid 
reference system the effects of the important nonlocal 
non-rigid electronic polarization processes in terms of lo- 
calized CF's and DF's have been studied. The latter lo- 
calized screening effects and not a homogenous electron 
gas screening allows for a quantitative calculation of the 
phonon anomalies and the dispersion. By the present 
investigations our earlier studies have been extended to 
the underdoped and overdoped state. Predictions of the 
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FIG. 9: Phonon density of states for LaCuO as calculated with the model for the insulating state (a), the underdoped state 
(b), the optimally doped state (c) and the overdoped state (d), respectively. 



dispersion in the latter case have been made. 

An important aspect of the calculations, besides the 
efforts to explore the phonon dynamics and the EPI in 
the HTSC's as strongly correlated systems, where these 
topics presently are not well understood, in particular 
in the insulating- and underdoped state, is to give addi- 
tional support to our modeling of the electronic state. 
This is achieved in terms of consecutive orbital selec- 
tive incompressibility-compressibility transitions. Such 
a modeling of the normal state qualitatively also pro- 
vides a possible picture of the superconducting state that 
points from a phase ordering transition in the under- 
doped regime to a BCS-like pairing transition for higher 



doping. Optimal Tc should mark a gradual crossover 
point. The localization-delocalization crossover in our 
approach corresponding to a phase ordering-pair order- 
ing crossover in the superconducting state is essentially 
related to the increase of an occupation with doping of 
the delocalized Cu4s orbitals progressively mixing into 
the Cu3d-02p dominated electronic state of the under- 
doped material. Simultaneously an emerging quasiparti- 
cle picture with a well defined Fermi surface is consistent 
with our modeling upon doping. Finally, the strong and 
material specific nonlocal EPI effects found in the cal- 
culations demonstrate the important role of the phonons 
for the physics of the cuprates. 
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